FTUV/09-1111 
IFIC/09-53 
LPT-Orsay-09-90 



Hadron structure in r —> KKttv t decays 



D. Gomez Dumm 1 , P. Roig 2 , A. Pich 3 , J. Portoles 3 



x ) IFLP, CONICET - Dpto. de Ffsica, Universidad Nacional de La Plata, 
C.C. 67, 1900 La Plata, Argentina 

2 ) Laboratoire de Physique Theorique (UMR 8627), Universite de Paris-Sud XI, 
Batiment 210, 91405 Orsay cedex, France 

3 ) Departament de Ffsica Teorica, IFIC, CSIC — Universitat de Valencia, 
Edifici d'Instituts de Paterna, Apt. Correus 22085, E-46071 Valencia, Spain 



Abstract 

We analyse the hadronization structure of both vector and axial- vector currents leading 
to r — > KKirv T decays. At leading order in the 1/Nc expansion, and considering only 
the contribution of the lightest resonances, we work out, within the framework of the res- 
onance chiral Lagrangian, the structure of the local vertices involved in those processes. 
The couplings in the resonance theory are constrained by imposing the asymptotic be- 
haviour of vector and axial-vector spectral functions ruled by QCD. In this way we pre- 
dict the hadron spectra and conclude that, contrarily to previous assertions, the vector 
contribution dominates by far over the axial-vector one in all KKir charge channels. 

PACS : 11.15.Pg, 12.38.-t, 12.39.Fe 
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1 Introduction 



Hadron decays of the tau lepton provide a prime scenario to study the hadronization 
of QCD currents in an energy region settled by many resonances. This task has a twofold 
significance. First, the study of branching fractions and spectra of those decays is a major goal 
of the asymmetric B factories (BABAR, BELLE). These are supplying an enormous amount of 
quality data owing to their large statistics, and the same is planned for the near future at tau- 
charm factories such as BES-III. Second, the required hadronization procedures involve QCD 
in a non-perturbative energy region (E < M T ~ 1.8 GeV) and, consequently, these processes 
are a clean benchmark, not spoiled by an initial hadron state, where we can learn about the 
treatment of strong interactions when driven by resonances. 

Analyses of tau decay data involve matrix elements that convey the hadronization of the 
vector and axial-vector currents. At present there is no determination from first principles 
of those matrix elements as they involve strong interaction effects in its non-perturbative 
regime. Therefore we have to rely in models that parameterize the form factors that arise 
from the hadronization. A relevant one is the so-called Kuhn-Santamaria model (KS) [T] 
that, essentially, relies on the construction of form factors in terms of Breit-Wigner functions 
weighted by unknown parameters that are extracted from phenomenological analyses of data. 
This procedure, that has proven to be successful in the description of the irnn final state, has 
been employed in the study of many two- and three-hadron tau decays [2HS]- The ambiguity 
related with the choice of Breit-Wigner functions [TJIB] is currently being exploited to estimate 
the errors in the determination of the free parameters. The measurement of the KKir spectrum 
by the CLEO Collaboration [7] has shown that the parameterization described by the KS model 
does not recall appropriately the experimental features keeping, at the same time, a consistency 
with the underlying strong interaction theory [8]. The solution provided by CLEO based in 
the introduction of new parameters spoils the normalization of the Wess-Zumino anomaly, 
i.e. a specific prediction of QCD. Indeed, arbitrary parameterizations are of little help in the 
procedure of obtaining information about non-perturbative QCD. They may fit the data but 
do not provide us hints on the hadronization procedures. The key point in order to uncover 
the inner structure of hadronization is to guide the construction of the relevant form factors 
with the use of known properties of QCD. 

The TAUOLA library j9] is, at present, a key tool that handles analyses of tau decay data. 
Though originally it comprehended assorted versions of the KS model only, it has been opened 
to the introduction of matrix elements obtained with other models. Hence it has become an 
excellent tool where theoretical models confront experimental data. This or analogous libraries 
are appropriate benchmarks where to apply the results of our research. 

At very low energies (E <C M p , being M p the mass of the p(770) resonance) the chiral 
symmetry of massless QCD rules the construction of an effective field theory that allows a 
perturbative expansion in momenta (p) and light quark masses (m), as (p 2 , M^)/A 2 , being 
A x ~ AtcF ~ M p the scale that breaks the chiral symmetry; here M n is the pion mass and F 
is the decay constant of the pion. Indeed Chiral Perturbation Theory (%PT) [TU] drives the 
hadronization of QCD currents into the lightest multiplet of pseudoscalar mesons, n, K and r\. 
The application of this framework to the study of pion decays of the tau lepton was carried out 
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in Ref. [TT] though, obviously, it can only describe a tiny region of the available phase space. 
It is clear that, whatever the structure given to the form factors in the region of resonances, 
it should match the chiral constraint in its energy domain. In fact the parameterizations with 
Breit-Wigner functions that are at the center of the KS model fail to fulfill that condition 
already at 0(p A ) in the chiral expansion [T2"|[l"3]. 

Our knowledge of QCD in the energy region of the resonances is pretty poor. Contrarily 
to the very low energy domain, we do not know how to construct a dual effective field theory 
of strong interactions for E ~ M p . There is a tool, though, that could shed light on the 
appropriate structure of a Lagrangian theory that we could use. This is yielded by the large- 
Nc limit of SU(Nc) QCD [14] , which introduces an expansion in inverse powers of the number 
of colours Nc- The essential idea relevant for our goal that comes out from that setting is 
that at leading order in the expansion, i.e. Nc —> oo, any amplitude is given by the tree 
level diagrams generated by a local Lagrangian with an spectrum of infinite zero-width states. 
This frame, as we will see, can be used to establish a starting point in the study of the hadron 
resonance region and, consequently, in the hadron decays of the tau lepton. The setting recalls 
the role of the resonance chiral theory [T5JII6] that can be better understood in the light of 
the large- N c limit [IMS]. 

At high energies (E ^> M p ), where the light-flavoured continuum is reached, perturbative 
QCD is the appropriate framework to deal with the description of strong interaction of partons. 
In particular, a well known feature of form factors of QCD currents is their smooth behaviour 
at high transfer of momenta [20], thus it is reasonable to expect that the form factors match 
this behaviour above the energy region of the resonances. Another related tool is the study of 
the Operator Product Expansion (OPE) of Green functions of QCD currents that are order 
parameters of the chiral symmetry breaking. It is possible to evaluate these Green functions 
within a resonance theory, and then perform a matching with their leading term of the OPE 
expansion at high transfers of momenta [2TH2T] . In general, the information coming from high 
energies is important to settle a resonance Lagrangian. It is reasonable to assume that the 
effective couplings collect information coming from the energy region above the resonances, 
hence the described procedure should help to determine the corresponding coupling constants. 
Indeed, this approach has proven to be capable of that task [27]. 

In Ref. [13] we considered all mentioned steps in order to analyse the nirir hadron final 
state in the decay of the tau lepton. Here we continue that undertaking by considering the 
KKtt channels that, as mentioned above, do not fit well within the KS model and the present 
TAUOLA setup. Contrarily to the ttutt final state, which is dominated by the hadronization 
of the axial-vector current, r — > KKttu t decays receive contributions from both vector and 
axial-vector currents. Indeed, one of the goals of our work is to find out the relative weight 
of those contributions. Fortunately we will be assisted in this task by the recent analysis of 
e + e~ — > KKn cross-section by BABAR |28j where a separation between isoscalar and isovector 
channels has been performed. Hence we will be able to connect both processes through CVC. 

In Section 2 we introduce the observables to be considered and the framework settled by 
the procedure sketched above. Then the amplitudes for KK7r decay channels are evaluated in 
Section 3. An analysis of how we can get information on the resonance couplings appearing 
in the hadronization of the currents is performed in Section 4. Finally we explain our results 
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in Section 5, and our conclusions are pointed out in Section 6. Four technical appendices 
complete our exposition. 



2 Theoretical framework 



The hadronization of the currents that rule semileptonic tau decays is driven by non- 
perturbative QCD. As mentioned in the Introduction, our methodology stands on the con- 
struction of an action, with the relevant degrees of freedom, led by the chiral symmetry and 
the known asymptotic behaviour of form factors and Green functions driven by large Nc QCD. 
We will limit ourselves to those pieces of the action that are relevant for the study of decays 
of the tau lepton into three pseudoscalar mesons. Hence we will need to include both even- 
and odd-intrinsic parity sectors. 

The large Nc expansion of SU(Nc) QCD implies that, in the Nc — > oo limit, the study 
of Green functions of QCD currents can be carried out through the tree level diagrams of 
a Lagrangian theory that includes an infinite spectrum of non-decaying states [H]. Hence 
the study of the resonance energy region can be performed by constructing such a Lagrangian 
theory. The problem is that we do not know how to implement an infinite spectrum in a model- 
independent way. However, it is well known from the phenomenology that the main role is 
always played by the lightest resonances. Accordingly it was suggested in Refs. [15j[T6] that 
one can construct a suitable effective Lagrangian involving the lightest nonets of resonances 
and the octet of Goldstone bosons states (tt, K and 77). This is indeed an appropriate tool to 
handle the hadron decays of the tau lepton. The guiding principle in the construction of such 
a Lagrangian is chiral symmetry. When resonances are integrated out from the theory, i.e. 
one tries to describe the energy region below such states (E <C M p ), the remaining setting is 
that of xPT, to which now we turn. 

The very low-energy strong interaction in the light quark sector is known to be ruled by the 
SU(3)l <8> SU(3)r chiral symmetry of massless QCD implemented in %PT. The leading even- 
intrinsic-parity 0(p 2 ) Lagrangian, which carries the information of the spontaneous symmetry 
breaking of the theory, is : 

F 2 



,(2) 
- X PT 



■(v"+x+> 



(1) 



where 



= i[u ] (d p - ir^)u - u(d^ - ii^v)] 



x± = u'xu'±ux'u , x = 2B (s + ip) , (2) 
and (...) is short for a trace in the flavour space. The Goldstone octet of pseudoscalar fields 
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is realized non-linearly into the unitary matrix in the flavour space 

u(tp) = expj-^$(x)J , (4) 

which under chiral rotations transforms as 

u(<p) ->■ 9ru{v) h(g,(pY = h(g,ip)u(ip)gl , (5) 

with g = (gL,gR) £ <S77(3)l <8> S77(3)r and h(g,<p) G 577(3) y. External hermitian matrix 
fields r M , s and p promote the global S{7(3)l ® S77(3)r symmetry to a local one. Thus, 
interactions with electroweak bosons can be accommodated through the vector t> M = (r fl +£ li )/2 
and axial-vector a M = (r^—£ fJi )/2 fields. The scalar field s incorporates explicit chiral symmetry 
breaking through the quark masses taking s = M. + . . ., with M. = diag(m u , md, m s ) and, 
finally, at lowest order in the chiral expansion F = F % = 92.4 MeV is the pion decay constant 
and B F 2 = -(0|^|0) . 

The leading action in the odd-intrinsic-parity sector arises at 0(p A ). This is given by the 
chiral anomaly [25] and explicitly stated by the Wess-Zumino-Witten Z wz [v,a] functional 
that can be read in Ref. (30]. This contains all anomalous contributions to electromagnetic 
and semileptonic meson decays. 

It is well known [11)|[1"9"] that higher orders in the chiral expansion, i.e. even-intrinsic-parity 
CyPT with n > 2, bring in the information of heavier degrees of freedom that have been 
integrated out, for instance resonance states. As our next step intends to include the latter 
explicitly, to avoid double counting issues we will not consider higher orders in x?T. As we 
comment below, in order to fulfill this procedure — at least, up to 0(p 4 ) — it is convenient 
to use the antisymmetric tensor representation for the J = 1 fields. Analogously, additional 
odd-intrinsic-parity amplitudes arise at 0(p 6 ) in xPT, either from one-loop diagrams using 
one vertex from the Wess-Zumino-Witten action or from tree-level operators |31j . However we 
will assume that the latter are fully generated by resonance contributions [23] and, therefore, 
will not be included in the following. 

The formulation of a Lagrangian theory that includes both the octet of Goldstone mesons 
and £7(3) nonets of resonances is carried out through the construction of a phenomenologi- 
cal Lagrangian [32] where chiral symmetry determines the structure of the operators. Given 
the vector character of the Standard Model (SM) couplings of the hadron matrix elements 
in t decays, form factors for these processes are ruled by vector and axial-vector resonances. 
Notwithstanding those form factors are given, in the r — > PPPu T decays, by a four-point Green 
function where other quantum numbers might play a role, namely scalar and pseudoscalar res- 
onances [33]. However their contribution should be minor for r — > KK-nv T . Indeed the lightest 
scalai0, namely / (980), couples dominantly to two pions, and therefore its role in the KKn 
final state should be negligible. Heavier flavoured or unflavoured scalars and pseudoscalars 
are at least suppressed by their masses, being heavier than the axial- vector meson ai(1260) 
(like J\"q(1430) that couples to Kir). In addition the couplings of unflavoured states to KK 

1 As we assume the Nc — > oo limit, the nonet of scalars corresponding to the /o(600) is not considered. This 
multiplet is generated by rescattering of the ligthest pseudoscalars and then subleading in the 1/Nc expansion. 
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(scalars) and KKtt (pseudoscalars) seem to be very small [35]. Thus in our description we 
include J = 1 resonances only 0, and this is done by considering a nonet of fields [15] : 

1 9 

R = —7= J2 K <pR,i , (6) 
V 2 i= i 

where R = V, A, stands for vector and axial- vector resonance states. Under the SU (3) l ® 
SU(3)r chiral group, R transforms as : 

R -> h{g,<p)Rh(jg,(pf . (7) 

The flavour structure of the resonances is analogous to that of the Goldstone bosons in Eq. (|3]). 
We also introduce the covariant derivative 

V^X = d„X + [T„,X\ , (8) 
r M = - [u\df, - ir^u + u(d^ - ii^u 1 ] , 

acting on any object X that transforms as R in Eq. ([7]), like w M and x±- The kinetic terms 
for the spin 1 resonances in the Lagrangian read : 

i /if 2 
£kin = -~{V X R X ,V V R^ ) + -f ( R^BT ) , R = V,A, (9) 

My, Ma being the masses of the nonets of vector and axial- vector resonances in the chiral 
and large-./Vc' limits, respectively. Notice that we describe the resonance fields through the 
antisymmetric tensor representation. With this description one is able to collect, upon in- 
tegration of resonances, the bulk of the low-energy couplings at 0{j> A ) in without the 
inclusion of additional local terms [27J. In fact it is necessary to use this representation if 
one does not include the C^ PT in the Lagrangian theory. Though analogous studies at higher 

chiral orders have not been carried out, we will assume that no £^pt w ^ n ti = 4,6, ... in the 
even-intrinsic-parity and n — 6,8, ... in the odd-intrinsic-parity sectors need to be included in 
the theory. 

The construction of the interaction terms involving resonance and Goldstone fields is driven 
by chiral and discrete symmetries with a generic structure given by : 

O t ~ (R 1 R 2 ...R jX {n) (<p)), (10) 

where X^iff) * s a chiral tensor that includes only Goldstone and auxiliary fields. It transforms 
like R in Eq. (J7J) and has chiral counting n in the frame of x?T. This counting is relevant 
in the setting of the theory because, though the resonance theory itself has no perturbative 
expansion, higher values of n may originate violations of the proper asymptotic behaviour of 
form factors or Green functions. As a guide we will include at least those operators that, 

2 If the study of these processes requires a more accurate description, additional resonances could also be 
included in our scheme. 
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contributing to our processes, are leading when integrating out the resonances. In addition 
we do not include operators with higher-order chiral tensors, x (v 9 ); that would violate the 
QCD asymptotic behaviour unless their couplings are severely fine tuned to ensure the needed 
cancellations of large momenta. In the odd-intrinsic-parity sector, that gives the vector form 
factor, this amounts to include all (Rx^) and (RRx^) terms. In the even-intrinsic-parity 
couplings, giving the axial- vector form factors, these are the terms (Rx^)- However previous 
analyses of the axial- vector contributions [131123] show the relevant role of the (RRx^) terms 
that, accordingly, are also considered hereQ 

We also assume exact SU(3) symmetry in the construction of the interacting terms, i.e. 
at level of couplings. Deviations from exact symmetry in hadronic tau decays have been 
considered in Ref. [31]. However we do not include SU(3) breaking couplings because we are 
neither considering next-to-leading corrections in the 1/Nc expansion. 

The lowest order interaction operators linear in the resonance fields have the structure 
(Rx^('-p))- There are no odd-intrinsic-parity terms of this form. The even-intrinsic-parity 
Lagrangian includes three coupling constants [U] : 

£2 = ^OWH, (ii) 

where f± u = uF^u^ ± u^F^u and F£ U L are the field strength tensors associated with the 
external right- and left-handed auxiliary fields. All coupling parameters Fy, Gy and Fa are 
real. 

The leading odd-intrinsic-parity operators, linear in the resonance fields, have the form 
(Rx^\<-p)) ■ We will need those pieces that generate : i) the vertex with one vector resonance 
and three pseudoscalar fields; ii) the vertex with one vector resonance, a vector current and 
one pseudoscalar. The minimal Lagrangian with these features is : 

i=l v i=l v 

where Qi and c, are real adimensional couplings, and the operators read 
1/ VPPP terms 

^vppp = «W ( yflu (ru\-^)) , 
^vppp = ie^ap {V^ (u y h orf u f) -u f) h or 'u 7 )) , 
Ovppp = W ({V^,u a u^} x-) , 



3 Operators {R\^) that are non-leading and have a worse high-energy behaviour, are not included in the 
even-intrinsic-parity contributions as they have not played any role in previous related analyses. 
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5 YPPP = e^g (u a V^u p X -) > (13) 
with ^ = V m m,, + V u u^ and 
2/ VJP terms [23] 

^ JP = w({^,/r}v^ ff ) , 

3 VJP = ie, upa ({V^JF} X -} , 

o YJP = £ ^({v tt y^/fK) , 

Owp = w({v^,/fK) • (14) 

Notice that we do not include analogous pieces with an axial-vector resonance, that would 
contribute to the hadronization of the axial- vector current. This has been thoroughly studied 
in Ref. [13] (see also Ref. [37]) in the description of the r — > -n-n-nv T process and it is shown that 
no (Ax^(<p)) operators are needed to describe its hadronization. Therefore those operators 
are not included in our minimal description of the relevant form factors. 

In order to study tau decay processes with three pseudoscalar mesons in the final state one 
also has to consider non-linear terms in the resonance fields. Indeed the hadron final state in 
r — > PPPu T decays can be driven by vertices involving two resonances and a pseudoscalar 
meson. The structure of the operators that give those vertices is (RiR^x^ '(¥>)) > an d has been 
worked out before [13j|2l]. They include both even- and odd- intrinsic-parity terms : 

5 4 

CT = X * °VAP + E ^ ^VVP , (15) 

i=l i=l 

where Aj, and di are unknown real adimensional couplings. The operators O l KKP are given by : 
1 / VAP terms 

Ol AP = ([V, A^x-) , 

2 VAP = i{[V^,A„ a )h«) , (16) 

<?vap = i([V»V» v ,A™]u a ) , 

Ol AP = i([V a V^,A^]u a ) . 
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2/ VVP terms 

0vvp = Wii^.^lV) , 
Ovvp = i£ M ,p CT ({^,^}X-) > 

^vvp = w({v^,^ Q K) • (17) 

We emphasize that C 2 R is a complete basis for constructing vertices with only one pseudoscalar 
meson; for a larger number of pseudoscalars additional operators might be added. As we are 
only interested in tree-level diagrams, the equation of motion arising from C x p T in Eq. (JT|) has 
been used in and C 2 R to eliminate superfluous operators. 

Hence our theory is given by the Lagrangian : 

£r x t = C>xpt + ^km + £2 + £2 + + ^2 R - ( 18 ) 

It is important to point out that the resonance theory constructed above is not a theory of 
QCD for arbitrary values of the couplings in the interaction terms. As we will see later on, 
these constants can be constrained by imposing well accepted dynamical properties of the 
underlying theory. 

3 Vector and axial- vector currents in r — > KKtt v t decays 

The decay amplitude for the r — > KKkv t decays can be written in the Standard Model as 

M = - -| V ud u VT Y (1 - 7s) u T T M , (19) 

where V u d is an element of the Cabibbo-Kobayashi-Maskawa matrix and T M is the hadron 
matrix element of the participating — QCD current : 

Tfi = (Kip,) K(p 2 ) 7r(p 3 ) I (V, - Afj) e iC ^ | 0) . (20) 

The hadron tensor can be written in terms of four form factors F,, F 2 , F 3 and F 4 as [M] : 

T M = VfFt + V£F 2 + V£F 3 + Q^F 4 , (21) 

where = + p 2 + p 3 and 

Vlp = y9»»-~^j (pa-pi)" > Vin = ( K 9^~~^qij (ps-pi^j 

V 3 » = ie^plplpl . (22) 
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There are three different charge channels for the KKn decays of the r lepton, namely 
K + {p + ) K~{pJ) 7T~(p n ), K°(po)K (p ) 7r~(j»7r) and K~(p^) K°(p ) n°(p n ). The definitions of 
Eq. fl22|) correspond to the choice p 3 = p n in all cases, and : (px,P2) = (p+,P-) f° r the K + K~ 
case, ipiiP-i) = (Po>Po) for K° K and (^1,^2) = (p_,po) for K~ K°. In general, form factors 
Fi are functions of the kinematical invariants : Q 2 , s = (pi +P2) 2 and t = (pi +P3) 2 . F x and F 2 
originate from the axial- vector current, while F 3 follows from the vector current. All of them 
correspond to spin-1 transitions. The F 4 pseudoscalar form factor stems from the axial-vector 
current, and corresponds to a spin-0 transition. It is seen that this form factor vanishes in 
the chiral limit, therefore its contribution is expected to be heavily suppressed, and both the 
spectrum and the branching ratio of tau decays into three pseudoscalar mesons is dominated 
by J = 1 transitions, especially in the Cabibbo-allowed modes. 

The Q 2 -spectrum is given by : 



dT_ _ G 2 F \V ud \ 2 (K_X , , , 
dQ 2 ~ 128(2tt)5M t \Q 2 ' ' (l " 



Wsa + Ui + 2^)(W a + W b ) 



(23) 



where the hadron structure functions, introduced in Ref. [36J, are : 
W A = -(Vl l F 1 + VfF 2 )(V lll F 1 + V 2fl F 2 )* , 

W B = \[stu + (m 2 K - ml) (Q 2 - m 2 K ) s + m 2 K (2m 2 n - Q 2 ) Q 2 - m 2 K ml] \F 3 \ 2 , 

W SA = (Q^F^Q.F:) = Q 2 \F A \ 2 , (24) 

where u = Q 2 — s — 1 + 2m 2 K + m 2 . The phase-space integrals extend over the region spanned 
by the hadron system with a center-of-mass energy \[Q 2 : 

dsdt = ds dt, (25) 

J Arn? K Jt-{s) 

with 

t±(s) = -L{(Q 2 -miy - [\V 2 (Q 2 ,s,ml) T A 1/2 {m 2 K ,m 2 K} s)f} , (26) 

and X(a,b,c) = (a + b — c) 2 — Aab. We have neglected here the v T mass, and exact isospin 
symmetry has been assumed. 

The general structure of the form factors, within our model, arises from the diagrams 
displayed in Fig. [TJ This provides the following decomposition : 

F = Ff + Ft + F™, z = l,..., (27) 

where Ff is given by the Lagrangian [topologies a) and b) in Fig. [I], and the rest are 
the contributions of one [Fig. lc), d) and e)] or two resonances [Fig. 1/)]. 
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a) 



b) 



c) 



d) 









/) 



Figure 1: Topologies contributing to the final hadron state in r — >■ KKir v T decays in the 
Nq — > oo limit. A crossed circle indicates the QCD vector or axial-vector current insertion. A 
single line represents a pseudoscalar meson (K, tc) while a double line stands for a resonance 
intermediate state. Topologies b) and e) only contribute to the axial-vector driven form factors, 
while diagram d) arises only (as explained in the text) from the vector current. 



3.1 Form factors in r 



K+K 7T" v T and r if° if tt" 



In the isospin limit, form factors for the r — > K + K it v t and r — >■ K°K°tt v t decays 
are identical. The explicit expressions for these are : 

y/2 



px 



3F ' 
6 F 3 
2 FaGv Q 2 



A R (Q 2 ,s,u,m 2 K , ml,m 2 K ) B n (s,u,m 2 K ,m 2 w ] 



M 2 -s 



M 2 K ,-t 



(28) 



3 F 3 Ml - Q 2 



A RR (Q 2 , s, u, m 2 K} ml,m 2 K ) 



M 2 - s 



+ 



B RR (Q 2 , s, u, t, m 2 K , ml, m 



Ki 



M 2 K * -t 



where the functions A n , B R , A RR and B RR are defined in Appendix |A] The dependence of the 
form factors with t follows from the relation u = Q 2 — s — 1+ 2m 2 K + m 2 . Moreover resonance 
masses correspond to the lowest states, M p = M p ( 770 ), M K * = M K *^ 2 ) and M ai = M ai ( 12 6o)0- 
Analogously the F 2 form factor is given by : 



px 



F x 



4 Resonance masses and widths within our approach are discussed in subsect. 3.3. 
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pa/ + \ FyGy 



Fr(s,t) 



6 F 3 
2 i^Gy Q 2 



B B (t,u,m 2 K , m 2 K ) A K (Q 2 ,t,u,m 2 K ,m 2 K ,ml) 



M 2 - s 



M 2 -t 



3 F 3 M 2 -Q 2 



B KR (Q 2 , t, u, s, m 2 K , m 2 K , m 



(29) 



M 2 -s 



+ 



A RR (Q 2 , t, u, m 2 K , m 2 K , ml) 



M%.-t 



The F 3 form factor arises from the chiral anomaly and the non-anomalous odd-intrinsic-parity 
amplitude. We obtain : 



r 3 



F*(s,t) 



N C V2 
12n 2 F 3 

AGy 
My F 3 



C R (Q 2 ,s. nr K . nr K .niz ) \ sur t) v 

+ COS 2 6y 



2 1 — y / 2tan0y 



M 2 - s 

2F V D*(Q 2 ,s,ty 
G v M 2 -Q 2 . 



+ 



1 + v^COt By 

M 2 -s 

C R (Q 2 ,t,m 2 K ,ml,m 2 K ) 
M 2 K *-t 



(30) 



Fr(s,t) 



AV2 Fv Gv ' 



F 3 M 2 - Q 2 



C RR (Q 2 , S ,m 2 ) sin 2 0, 



1 + v^cot 9\ 
M 2 - s 



+ cos 2 01 



1 - V2tan0 v 
M 2 -s 



C™(Q 2 ,t,m 2 K ) 
M%.-t 



where C R , D K and C RR are defined in Appendix [A^ and 9y is the mixing angle between the 
octet and singlet vector states u$ and cu that defines the mass eigenstates cu(782) and 0(1020) : 



U) 



cos 9y — sin 9y 
sin 9y cos 9y 



UJ$ 

CJ 



(31) 



For numerical evaluations we will assume ideal mixing, i.e. 9 V = tan-^l/v 7 ^). In this case 
the contribution of the 0(1020) meson to F3 vanishes. 

Finally, though we have not dwelled on specific contributions to the F4 form factor, we 
quote for completeness the result obtained from our Lagrangian. Its structure is driven by the 
pion pole : 



F, 

F*(s,t) 



Ff + Ff, 



mz 



V2Fml - Q 1 



1 + 



nix — u 



Q 2 
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Gl7 



ml 



y/2 F 3 Q 2 (ml-Qi 



s(t — u) t(s 



u 



) - {m\ - ml)(Q 2 - m 2 K ) 



M 2 p 



Ml, 



• (32) 



3.2 Form factors in r — > K K° 7r° za 



The diagrams contributing to the r~ — > K~ K° tc° v t decay amplitude are also those in Fig. (TJ 
hence once again we can write Fi = Ff + F R + F RR + . . . . However, the structure of the 
form factors for this process does not show the symmetry observed in r — > KKttu t . We find : 

1 



if 
if(s,t) 



1 FyGy 



6 F 3 



B R (s,u,m 2 K ,ml) ^ A R (Q 2 , s, u, m 2 K , m 2 , m 2 K ) 



M 2 K ,-t 



M 2 p - s 



+ 



A n (Q 2 ,u,s,ml,m 2 K ,m 2 K ) 



V2F A G V Q 2 
3 F 3 M 2 -Q< 



M 2 K , - u 
B nn {Q 2 ,s,u,t,m 2 K ,ml,m 2 K ) 



M 2 K ,-t 



+ 2 



A RR (Q 2 , s, u, m 2 K , m 2 , m 2 K ) 



M 2 p - s 



+ 



A RR (Q 2 1 u 1 s,m 2 7r ,m 2 K ,m 2 K ) 



M 2 K , - u 



(33) 



Ff(8,t) 



1 FyGy 



6 F 3 



A R (Q 2 ,t,u,m 2 K ,m 2 K} m%)_ g B R {t,u } m 2 K} m 2 K ) 



M 2 K ,-t 



M 2 p - s 



A K {Q 2 ,u,t,m 2 K ,m 2 K ,ml) 



u 



F™(s,t) 



V2F A G V Q 2 
3 F 3 M 2 -Q> 



A RR (Q , t, u, m K} m K} m^) 



M 2 K ,-t 



+ 2 



B RIi (Q 2 ,t,u, s,m 2 K ,m 2 K ,ml 
M 2 p - s 
A RR {Q 2 ,u,t,m 2 K ,m 2 K ,m 2 n ) 
M 2 K ,-u 



The form factor driven by the vector current is given by : 
if = 



(34) 



F R (s,t) 



2V2G V 

My F 3 



C R (Q 2 ,t,m 2 K ,mlm 2 K ) C R (Q 2 ,u,m 2 K ,ml,m^) 2F V E R (t,u 



M 2 K *-t 



u 



Gy M 2 -Q< 



12 



F™(s,t) 



-4 



FyGy 



1 



f 3 m 2 - g 2 



C RR (Q 2 ,t,m^) C RR (g 2 ,n,m^) 



M 2 K ,-t 



M 2 K ,-u 



(35) 



with i? R defined in Appendix lAl 

Finally for the pseudoscalar form factor we have 

1 m\{t- u) 



F*(s,t) 



2FQ 2 (m 2 n -Q 2 ) ' 

ml \ t(s - u) - {m\ - ml)(Q 2 - m\) 2s(t-u) 



IG 2 V 



2 F 3 Q 2 {ml-Q 2 ) 



M%, -t 



+ 



M 2 - s 



u(s - 1) - (m 2 K - m\){Q 2 - m 2 K ) 
M£. -u 



(36) 



3.3 Features of the form factors 

Several remarks are needed in order to understand our previous results for the form factors 
related with the vector and axial-vector QCD currents analysed above : 

1 / Our evaluation corresponds to the tree level diagrams in Fig. [T] that arise from the 
N c —> oo limit of QCD. Hence the masses of the resonances would be reduced to 



Mi 



Mk* = Ms and Ma = M ai as they appear in the resonance 



v - 

Lagrangian (Q, i.e. the masses of the nonet of vector and axial- vector resonances in the 
chiral and large- Nc limit. However it is easy to introduce NLO corrections in the 1/Nc 
and chiral expansions on the masses by including the physical ones : M p , M^*, M^, M^ 
and M ai for the p(770), K*(892), w(782), 0(1020) and ai(1260) states, respectively, as 
we have done in the expressions of the form factors. In this setting resonances also have 
zero width, which represents a drawback if we intend to analyse the phenomenology of 
the processes : Due to the high mass of the tau lepton, resonances do indeed resonate 
producing divergences if their width is ignored. Hence we will include energy- dependent 
widths for the p(770), ai(1260) and if* (892) resonances, that are rather wide, and a 
constant width for the cj(782). This issue is discussed in Ref. [37] . 

In summary, to account for the inclusion of NLO corrections we perform the substitu- 
tions : 

1 1 

(37) 



M 2 R ~ q 2 M 2 hys - q 2 - i M phys T phys (q 2 ) ' 

where R = V, A, and the subindex phys on the right hand side stands for the correspond- 
ing physical state depending on the relevant Feynman diagram. 

2/ If we compare our results with those of Ref. [3], evaluated within the KS model, we 
notice that the structure of our form factors is fairly different and much more intricate. 
This is due to the fact that the KS model, i.e. a model resulting from combinations of 
ad hoc products of Breit- Wigner functions, does not meet higher order chiral constraints 
enforced in our approach. 
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3/ As commented above the pseudoscalar form factors F4 vanishes in the chiral limit. Indeed 
the results of Eqs. ( 13211361) show that they are proportional to m 2 , which is tiny compared 
with any other scale in the amplitudes. Hence the contribution of F4 to the structure of 
the spectra is actually marginal. 

4 QCD constraints and determination of resonance cou- 
pling constants 

Our results for the form factors Fi depend on several combinations of the coupling constants 
in our Lagrangian Cr x t in Eq. (fl8"]) . most of which are in principle unknown parameters. Now, 
if our theory offers an adequate effective description of QCD at hadron energies, the underlying 
theory of the strong interactions should give information on those constants. Unfortunately 
the determination of the effective parameters from first principles is still an open problem in 
hadron physics. 

A fruitful procedure when working with resonance Lagrangians has been to assume that 
the resonance region, even when one does not include the full phenomenological spectrum, 
provides a bridge between the chiral and perturbative regimes [27]. The chiral constraints 
supply information on the structure of the interaction but do not provide any hint on the 
coupling constants of the Lagrangian. Indeed, as in any effective theory [39], the couplings 
encode information from high energy dynamics. Our procedure amounts to match the high 
energy behaviour of Green functions (or related form factors) evaluated within the resonance 
theory with the asymptotic results of perturbative QCD. This strategy has proven to be 
phenomenologically sound [2TM2~7] |4T)] . and it will be applied here in order to obtain information 
on the unknown couplings. 

Two-point Green functions of vector and axial- vector currents Uv,a{c[ 2 ) were studied within 
perturbative QCD in Ref . [H] , where it was shown that both spectral functions go to a constant 
value at infinite transfer of momenta : 

^mU v>A (q 2 ) (38) 

By local duality interpretation the imaginary part of the quark loop can be understood as the 
sum of infinite positive contributions of intermediate hadron states. Now, if the infinite sum 
is going to behave like a constant at q 2 — > 00, it is heuristically sound to expect that each one 
of the infinite contributions vanishes in that limit. This deduction stems from the fact that 
vector and axial- vector form factors should behave smoothly at high q 2 , a result previously 
put forward from parton dynamics in Ref. [20] . Accordingly in the Nq — > 00 limit this result 
applies to our form factors evaluated at tree level in our framework. 

Other hints involving short- distance dynamics may also be considered. The analyses of 
three-point Green functions of QCD currents have become a useful procedure to determine 
coupling constants of the intermediate energy (resonance) framework [2TM25] . The idea is to 
use those functions (order parameters of the chiral symmetry breaking), evaluate them within 
the resonance framework and match this result with the leading term in the Operator Product 
Expansion (OPE) of the Green function. 
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In the following we collect the information provided by these hints on our coupling con- 
stants, attaching always to the Nq — > oo case [18] (approximated with only one nonet of vector 
and axial-vector resonances) : 

i) By demanding that the two-pion vector form factor vanishes at high q 2 one obtains the 
condition Fy Gy = F 2 involving the couplings in Eq. (11 II) 



ii) The first Weinberg sum rule [43] leads to F v — F\ = F 2 , and the second Weinberg sum 
rule gives F v M v = F\ M% [E] . 

iii) The analysis of the VAP Green function [23] gives for the combinations of couplings 
defined in Eq. (I A . 2 [) the following results : 



A' 



A" 



Mi 



2^2 F A G V 2 V2 M v " 
2G V -F v Ml - 2M 2 



2V2 F A 2^f2M v M A 

4A = A' + A", (39) 

where, in the two first relations, the second equalities come from using relations i) and 
ii) above. Here My and M A are the masses appearing in the resonance Lagrangian 
(Q. Contrarily to what happens in the vector case where My is well approximated 
by the p(770) mass, in Ref. [26] it was obtained M A = 998(49) MeV, hence M A differs 
appreciably from the presently accepted value of M ai . It is worth to notice that the two 
first relations in Eq. ( 1391) can also be obtained from the requirement that the J = 1 axial 
spectral function in r — >■ ?>ttv t vanishes for large momentum transfer [T3] . 

iv) Both vector form factors contributing to the final states KKn~ and K~K°7i° in tau 
decays, when integrated over the available phase space, should also vanish at high Q 2 . 
Let us consider H^(s,t,Q 2 ) = T^T^*, where can be inferred from Eq. ( 12TI) . Then 
we define Uy(Q 2 ) by : 

J dU 3 H^(s,t,Q 2 ) = (Q 2 g, u - U V (Q 2 ) , (40) 



where 



/dn 3 = J d ^ d ^^6HQ-Pi-P2-Pz)Hs-(Q-P3) 2 )Ht-(Q-P2) 2 ) 

dsdt. (41) 



TV 2 



4Q 2 

Hence we find that 



UviQ2) = lW / dsdt ^^M>^)> ( 42 ) 
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where the limits of integration are those of Eq. ( 125] 126]) . should vanish at Q — > oo. 
This constraint determines several relations on the couplings that appear in the F 3 form 
factor, namely : 

ci - c 2 + c 5 = , (43) 

N c F V M V 

Cl _ C2 _ C5 + 2ce = -jjj^^, (44) 

4 = -Jk-&, (45) 

3 192 vr 2 F 2 v ; 

i?i + 2# 2 -<7 3 = 0, (46) 

N c My 

92 = — k (47) 



192^ 



7T 



If these conditions are satisfied, Uy(Q 2 ) vanishes at high transfer of momenta for both 
KKtt~ and K~K°7i° final states. We notice that the result in Eq. ( 143]) is in agreement 
with the corresponding relation in Ref. [24], while Eqs. (144]) and ( 145]) do not agree with 
the results in that work. In this regard we point out that the relations in Ref. [24J, 
though they satisfy the leading matching to the OPE expansion of the (VVP) Green 
function with the inclusion of one multiplet of vector mesons, do not reproduce the right 
asymptotic behaviour of related form factors. Indeed it has been shown [22,26j that two 
multiplets of vector resonances are needed to satisfy both constraints. Hence we will 
attach to our results above, which we consider more reliable 0. 

v) An analogous exercise to the one in iv) can be carried out for the axial-vector form 
factors Fi and F 2 . We have performed such an analysis and, using the relations in i) 
and ii) above, it gives us back the results provided in Eq. (139]) for A' and A". Hence both 
procedures give a consistent set of relations. 

After imposing the above constraints, let us analyse which coupling combinations appearing in 
our expressions for the form factors are still unknown. We intend to write all the information 
on the couplings in terms of F, My and Ma- From the relations involving Fy, Fa and Gy we 
obtain : 



n 


Ml 


F 2 


Ml -My 


n 


My 


F 2 


Ml- My 


r<2 

Lry 


, My 


F 2 " 





(48) 

Moreover we know that Fy and Gy have the same sign, and we will assume that it is also the 
sign of Fa- Together with the relations in Eq. ( 139]) this determines completely the axial-vector 



5 One of the form factors derived from the (VVP) Green function is J r 7r7 * 7 ((7 2 ), that does not vanish at 
high q 2 with the set of relations in Ref. [24]. With our conditions in Eqs. (|44I45[) the asymptotic constraint 
on the form factor can be satisfied if the large- Nc masses, Ma and My, fulfill the relation 2M\ — 3M V . It is 
interesting to notice the significant agreement with the numerical values for these masses mentioned above. 
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form factors F\^- Now from Eqs. (14311471) one can fix all the dominant pieces in the vector form 
factor F 3 , i.e. those pieces that involve factors of the kinematical variables s, t or Q 2 . The 
unknown terms, that carry factors of m 2 or m 2 K , are expected to be less relevant. They are 
given by the combinations of couplings : c\ + c 2 + 8 c 3 — c 5 , d\ + 8 d 2 , c 4 , g 4 and g 5 . However 
small they may be, we will not neglect these contributions, and we will proceed as follows. 
Results in Ref. [23] determine the first and the second coupling combinations. As commented 
above the constraints in that reference do not agree with those we have obtained by requiring 
that the vector form factor vanishes at high Q 2 . However, they provide us an estimate to 
evaluate terms that, we recall, are suppressed by pseudoscalar masses. In this way, from a 
phenomenological analysis of u — > 7r + 7r~7r° (see Appendix [B]) it is possible to determine the 
combination 2 g^+g^. Finally in order to evaluate c 4 and g 4 we will combine the recent analysis 
of cr (e + e~ — > KKtt) by BABAR [28] with the information from the r — > KK-nv T width. 



4.1 Determination of C4 and 



The separation of isoscalar and isovector components of the e + e~ — > KKtt amplitudes, carried 
out by BABAR |28j, provides us with an additional tool for the estimation of the coupling 
constant c 4 that appears in the hadronization of the vector current. Indeed, using SU{2)i 
symmetry alone one can relate the isovector contribution to a (e + e- K~K°ir + ) with the 
vector contribution to T (t~ — > K Q K~it v T ) through the relation : 



f(Q 2 ) <ti=i (e + e- -> K-K\ + ) 



(49) 



f, 



where f(Q 2 ) is given in Appendix In this Appendix we also discuss other relations similar 
to Eq. ( 14"9"|) that have been used in the literature and we point out the assumptions on which 
they rely. 

Hence we could use the isovector contribution to the cross-section for the process e + e~ — > 
KsK ± n Zf determined by BABAR and Eq. (14"9~|) to fit the c 4 coupling that is the only still unde- 
termined constant in that process. However we have to take into account that our description 
for the hadronization of the vector current in the tau decay channel does not, necessarily, 
provide an adequate description of the cross-section. Indeed the complete different kinematics 
of both observables suppresses the high-energy behaviour of the bounded tau decay spectrum, 
while this suppression does not occur in the cross-section. Accordingly, our description of the 
latter away from the energy threshold can be much poorer. As can be seen in Fig. [2] there is 
a clear structure in the experimental points of the cross-section that is not provided by our 
description. 

Taking into account the input parameters quoted in Eq. fID.lj) we obtain : c 4 = — 0.047 ± 
0.002. The fit has been carried out for the first 6 bins (up to E cm ~ 1.52 GeV). This result 
corresponds to x 2 /dof = 0.3 and the displayed error comes only from the fit. 

We take into consideration now the measured branching ratios for the KKtt channels of 
Tabled] in order to extract information both from c 4 and g 4 . We notice that it is not possible 
to reconcile a prediction of the branching ratios of r — > KKixv T and r — > K~ K°ir Q v T in spite 
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Figure 2: Comparison of the experimental data [28] with the theoretical prediction for the cross- 
section of the isovector component of e + e~ — > K*(892)K — > KsK^tt^ process, for different values 
of the C4 coupling. The x 2 values are associated to the first 6 data points only. 



of the noticeable size of the errors shown in the Table [TJ Considering that the second process 
was measured long ago and that the t~ — > K + K~iz~v T decay has been focused by both CLEO 
III and BABAR we intend to fit the branching ratio of the latter. For the parameter values : 

c 4 = -0.07 ±0.01, 

34 = -0.72 ±0.20, (50) 

we find a good agreement with the measured widths r(r~ — > K + K~ir~i/ T ) and T(r — > 
K~ K°Tt Q v T ) within errors (see Tabled]). Notice that the value of \c±\ is larger than that ob- 
tained from the fit to the e + e~ — > KsK^w^ data explained above. In Fig. [2] we show the first 
8 bins in the isovector component of e + e~ — > KsK ± it t and the theoretical curves for different 
values of the C4 coupling. As our preferred result we choose the larger value of C4 in Eq. ( 1501 . 
since it provides a better agreement with the present measurement of r(r~ — > K~ K°7t Q v T ). 
Actually, one can expect a large incertitude in the splitting of isospin amplitudes in the 
e + e~ — > KsK ± ir T cross-section (see Appendix [C]) . Taking into account this systematic error, 
it could be likely that the theoretical curve with C4 = —0.07 falls within the error bars for the 
first data points. 
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Source 


r(r" -> K + K-n-iy T ) 


r(r~ ->■ k k\-v t ) 


r(r" -> K-K°ir°is T ) 


rUb [ooj 
BABAR HI] 
CLEO III [7] 

Belle gS] 
Our prediction 


o.lUo (loo J 

3.049 (85) 

3.511 (245) 

3.465 (136) 
o 4+0.5 


O.400 (_ mUJ 
o 4+O.5 


o./OZ (ozlj 

o c+0.3 
z.j_ 2 



Table 1: Comparison of the measurements of partial widths (in units of 10 15 GeV) with our pre- 
dictions for the set of values in Eq. (|5U|) . For earlier references see |35j . 



5 Phenomenology of r — > KKnu T : Results and their 
analysis 

Asymmetric B-factories span an ambitious r programme that includes the determination of 
the hadron structure of semileptonic r decays such as the KKn channel. As commented in the 
Introduction the latest study of t~ — > K + K~tt~v t by the CLEO III Collaboration [7J showed 
a disagreement between the KS model, included in TAUOLA, and the data. Experiments with 
higher statistics such as BABAR and Belle should clarify the theoretical settings. 

For the numerics in this Section we use the values in Appendix [Dl At present no spectra 
for these channels is available and the determinations of the widths are collected in Table [TJ 
We also notice that there is a discrepancy between the BABAR measurement of r(r~ — > 
K + K~7T~u T ) and the results by CLEO and Belle. Within SU(2) isospin symmetry it is found 
that r(r~ — > K + K~tt~v t ) = r(r~ — > K°K ir~v T ), which is well reflected by the values 
in Table [1] within errors. Moreover, as commented above, the PDG data [35] indicate that 
r(r~ — > K~ i^°7r°z/ T ) should be similar to r(r~ — > KKixv T ). It would be important to obtain 
a more accurate determination of the — > K~ K°7T°b> T width (the measurements quoted by 
the PDG are rather old) in the near future. 

In our analyses we include the lightest resonances in both the vector and axial-vector 
channels, namely p(775), K*(892) and Gq(1260). It is clear that, as it happens in the r — > 
ixixixv T channel (see Ref. [37]), a much lesser role, though noticeable, can be played by higher 
excitations on the vector channel. As experimentally only the branching ratios are available 
for the KK7r channel we think that the refinement of including higher mass resonances should 
be taken into account in a later stage, when the experimental situation improves. 

In Figs. |3] and H] we show our predictions for the normalized M\ Klx — spectrum of the 
t~ — > K + K~-k~i> t and r _ — > K~ K°tt°v t decays, respectively. As discussed above we have 
taken C4 = — 0.07±0.01 and g± = —0.72 ±0.20 (notice that the second process does not depend 
on (74). We conclude that the vector contribution (r^) dominates over the axial- vector one 
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(Ta) in both channels 



rv 



= 0.16 ±0.05 , 

KKtt Ly 



r r ->• K+K-tt-v t ) 
= 0.18 ±0.04, r ^ — n V = 1-4 ±0.3 



(51) 

where the errors estimate the slight variation due to the range in c\ and g^. These ratios 
translate into a ratio of the vector current to all contributions of f v = 0.86 ± 0.04 for the 
KKtt channel and f v = 0.85 ± 0.03 for the K~K°tc one, to be compared with the result in 
Ref. [12], namely f v (KKir) = 0.20 ± 0.03. Our results for the relative contributions of vector 
and axial-vector currents deviate strongly from most of the previous estimates, as one can 
see in Table |2j Only Ref. [5] pointed already to vector current dominance in these channels, 
although enforcing just the leading chiral constraints and using experimental data at higher 
energies. We conclude that for all r — > KK-kv t channels the vector component dominates by 
far over the axial- vector one, though, as can be seen in the spectra in Figs. I3"f4| the axial-vector 
current is the dominant one in the very-low Q 2 regime. 

Next we contrast our spectrum for r~ — > K + K~7T~v T with that one arising from the KS 
model worked out in Refs. [31I1S]- This comparison is by no means straight because in these 
references a second and even a third multiplet of resonances are included in the analysis. As 
we consider that the spectrum is dominated by the first multiplet, in principle we could start 
by switching off heavier resonances. However we notice that, in the KS model, the p(1450) 
resonance plays a crucial role in the vector contribution to the spectrum. This feature depends 
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UVJ LL-L v < V... 


iv/r^ 


Our result 


6±2 


Kb model [oj 


0.6 — U. / 


KS model 06] 


0.4-0.6 


Breit-Wigner approach [5] 


~ 9 


CVC @2] 


0.20 ±0.03 


Data analysis [7J 


1.26 ±0.35 



Table 2: Comparison of the ratio of vector and axial- vector contribution for r — > KKirv T partial 
widths. The last two lines correspond to the r~ — > K + K~n~v T process only. Results in Ref. [16] 
are an update of Ref. [3J- The result of Ref. [32] is obtained by connecting the tau decay width with 
the CVC related e + e~ — > KsK^tt^ (see Appendix [C]). The analysis in [7J was performed with a 
parameterization that spoiled the chiral normalization of the form factors. 




Figure 4: Normalized Mj^^-spectra for r — > K K 7r v T . Notice the dominance of the axial- vector 
current at very low values of Q 2 . 



strongly on the value of the p(1450) width, which has been changed from Ref. [3] to Ref. jl6] [§. 
In Fig. [5] we compare our results for the vector and axial- vector contributions with those of 

6 Moreover within Ref. [3] the authors use two different set of values for the p(1450) mass and width, one 
of them in the axial- vector current and the other in the vector one. This appears to be somewhat misleading. 
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the KS model as specified in Ref. [46] (here we have switched off the seemingly unimportant 
K*(1A10)). As it can be seen there are large differences in the structure of both approaches. 
Noticeably there is a large shift in the peak of the vector spectrum owing to the inclusion 
of the p(1450) and p(1700) states in the KS model together with its strong interference with 
the p(770) resonance. In our scheme, including the lightest resonances only, the p(1450) and 
p(1700) information has to be encoded in the values of C4 and #4 couplings (that we have 
extracted in Subsection I4.ip and such an interference is not feasible. It will be a task for the 
experimental data to settle this issue. 
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Figure 5: Comparison between the normalized M^-^-^-spectra for the vector and axial- vector con- 
tributions to the t~ — > K^K^Ti^Vr channel in the KS model |46| and in our approach. 



In Fig. [6] we compare the normalized full M] (Kn spectrum for the r — > KKt^v t channels 
in the KS model [46J and in our scheme. The most important feature is the large effect of the 
vector contribution in our case compared with the leading role of the axial- vector part in the 
KS model, as can be seen in Fig. [5j This is the main reason for the differences between the 
shapes of M] (Kn spectra observed in Fig. [6j 



6 Conclusions 

Hadron decays of the tau lepton are an all-important tool in the study of the hadronization 
process of QCD currents, in a setting where resonances play the leading role. In particular 
the final states of three mesons are the simplest ones where one can test the interplay between 
different resonance states. At present there are three parameterizations implemented in the 
TAUOLA library to describe the hadronization process in tau decays. Two are based on 
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T -> K + K 7t V 
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Figure 6: Comparison between the normalized -spectra for r — > K + K ir v T in the KS 

model [36] and in our approach. 
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Figure 7: Comparison between the normalized M|-~ -spectra for r 
model [M] and in our approach. 
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experimental data. The other alternative, namely the KS model, though successfull in the 
account of the titttt final state, has proven to be unsuitable [7] when applied to the decays 
into KKtt hadron states. Our procedure, guided by large Nc, chiral symmetry and the 
asymptotic behaviour of the form factors driven by QCD, was already employed in the analysis 
of t — > 7T7T7rz/ r in Refs. |J3] and [37], which only concern the axial- vector current. Here we 
have applied our methodology to the analysis of the r — > KK-kv t channels where the vector 
current may also play a significant role. 

We have constructed the relevant Lagrangian involving the lightest multiplets of vector and 
axial-vector resonances. Then we have proceeded to the evaluation of the vector and axial- 
vector currents in the large-iVc limit of QCD, i.e. at tree level within our model. Though 
the widths of resonances are a next-to- leading effect in the 1/Nc counting, they have to be 
included into the scheme since the resonances do indeed resonate due to the high mass of the 
decaying tau lepton. We have been able to estimate the values of the relevant new parameters 
appearing in the Lagrangian with the exception of two, namely the couplings C4 and g±, which 
happen to be important in the description of r — > KKnv T decays. The range of values for 
these couplings has been determined from the measured widths r(r~ — > K + K~n~v T ) and 
r(r" -> K-K°tt°u t ). 

In this way we provide a prediction for the — still unmeasured — spectra of both processes. 
We conclude that the vector current contribution dominates over the axial-vector current, in 
fair disagreement with the corresponding conclusions from the KS model j3H] with which we 
have also compared our full spectra. On the other hand, our result is also at variance with 
the analysis in Ref. [32] • There are two all-important differences that come out from the 
comparison. First, while in the KS model the axial- vector contribution dominates the partial 
width and spectra, in our results the vector current is the one that rules both spectrum 
and width. Second, the KS model points out a strong interference between the p(770), the 
p(1450) and the p(1700) resonances that modifies strongly the peak and shape of the M KKn 
distribution depending crucially on the included spectra. Not having a second multiplet of 
vector resonances in our approach, we cannot provide this feature. It seems strange to us 
the overwhelming role of the p(1450) and p(1700) states but it is up to the experimental 
measurements to settle this issue. 

Even if our model provides a good deal of tools for the phenomenological analyses of ob- 
servables in tau lepton decays, it may seem that our approach is not able to carry the large 
amount of input present in the KS model, as the later includes easily many multiplets of 
resonances. In fact, this is not the case, since the Lagrangian can be systematically extended 
to include whatever spectra of particles are needed. If such an extension is carried out the 
determination of couplings could be cumbersome or just not feasible, but, on the same footing 
as the KS model, our approach would provide a parameterization to be fitted by the experi- 
mental data. The present stage, however, has its advantages. By including only one multiplet 
of resonances we have a setting where the procedure of hadronization is controlled from the 
theory. This is very satisfactory if our intention is to use these processes to learn about QCD 
and not only to fit the data to parameters whose relation with the underlying theory is unclear 
when not directly missing. 

We intend to follow our approach to analyse further relevant three pseudoscalar channels 
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along the lines explained in this article. 
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A Definitions in the expressions of form factors 

The results for the F\, F 2 and F 3 form factors in r — > KKttv t decays given in Eqs. ( )28|) . (129"]) . 
(|30|) . ( )33|) . (|34|) and (|35|) are expressed in terms of the following functions : 

A n (Q 2 ,x,y,ml,m 2 2 ,ml) = 3 x + m\ - m 2 3 + ( 1 - ^ j [2 Q 2 - 2 x - y + m 2 3 - ml] , 



B R (x,y,ml,ml) 




RR 



(Q 2 ,x,y,m\,ml,ml) 



(A' + A") (Sx + ml-ml) 




B RR (Q 2 ,x,y,z,m 2 1 ,m 2 2 ,m 2 3 ) 



2 (A' + A") (ml 



— m. 




C R (Q 2 ,x,ml,m 2 2 ,m 2 3 ) 



(ci - c 2 + c 5 ) Q 2 - (ci - c 2 - c 5 + 2c 6 ) x 
+ (ci + c 2 + 8c 3 - c 5 ) mj + 8 c 4 (ml - m\) , 

d 3 (Q 2 + x) + (di + 8d 2 - d 3 ) m 2 , 



C RR (Q 2 ,x,m 2 ) 



D R (Q 2 ,x,y) 



(gi + 2g2-g 3 ) (x + y)~2g 2 (Q 2 + m 2 K ) 

~(gi - 93) (3 m 2 K + ml) + 2 g 4 (m 2 K + ml) +2g 5 m 2 K , 



E R (x,y) 



(gi + 2g 2 -g 3 ) (x-y). 



(A.l) 
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Here u = Q 2 — s — t + m\ + m 2 , + m\ and F(x, y) = X" + A' x — Ao y, where Ao, A' and A" are 
combinations of the Aj couplings defined in Eq. ( ITS]) : 

-V2\ = 4A 1 + A 2 + ^ + A 5 , 

y/2X' = A 2 -A 3 + y + A 5 , (A.2) 

V2\" = A 2 -^-A 5 . 



B 2 g 4 + g 5 from — > 7r + 7r 7r° 

The process u — > ii + (ki) 7r~(/c 2 ) 7r°(/c3) provides us with an estimate for the combination of 
couplings 2 gi + g$. We will denote the polarization vector of the u as and use the kinematic 
invariants Sij = (ki + kj) 2 . 

The amplitude for this process has two contributions. The first one, mediated by the p(770) 
resonance was already studied in Ref. [21], where it was concluded that the contribution of a 
pure local amplitude was necessary to fulfill the phenomenological determinations. This piece 
can be obtained from our Lagrangian in Eq. ( IT3]) . The full result is given by : 

8 Gy 

m\{d x + 8d 2 - d 3 ) + (M 2 + s 12 )d 3 

TJl V 1^12 -> Si 3 | + |S 12 -> S 23 j 

where we have assumed ideal mixing between the states | ui 8 ) and | Ui) : 

|w> = y||wi) + y|(o; 8 ). (B.2) 

Using the experimental figure for BR(u; — > 7r + 7r~7r°) = 0.892 ± 0.007 (35], introducing the 
already known combinations of couplings as discussed in Sect. @]and taking Gy = Fy/F 2 with 
Fy given by Eq. (ID.lj) we find : 

2 # 4 + £ 5 = -0.60 ±0.02. (B.3) 
We will use this result to eliminate g$ in terms of g^, that remains unknown. 



C Relation between a (e + e — > KKtt) and r (r — > KKtiv t ) 

Using SU(2)i symmetry, one can derive several relations between exclusive isovector hadron 
modes produced in e + e~ collisions and those related with the vector current (F 3 form factor) 
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in t decays. In particular we find 



dQ : 



T (r - -> K°K-n°is T ) 



f(Q 2 ) cr I=1 (e+e" -> K-K\ + ) , 



where 



AV ' 128(2ir) 5 Af r \ Q* V 3 V M ? ) V » 2 I 



Analogously, one can also derive: 



F 3 



f(Q 2 ) [<ti=i (e+e" -> K+^tt") 

+ 2(T /=1 (e+e" -»■ K+K-Ti ) 



(C.l) 



(C.2) 



(C.3) 



Summing these equations one obtains : 



i-3 



/(Q 2 ) £ <x /=1 (e+e" 



(^TT)i) 



(C.4) 



where the sums run over all possible charge channels in each case. If isovector and isoscalar 
components were splitted for all channels Eq. ( 1C.4|) would allow us to fit the data using our 
vector form factors for r — > {KK^\v T . 

BABAR has managed to split the isoscalar and isovector components in the cross sections 
a (e+e" — > KsK ± ir T ) [28]. The 1 = 1 component of a (e+e~ — > KK tt) needs to be used, under 
the hypothesis of CVC, to obtain the spectral function of the processes r — > KKtxv t , and thus 
to help the extraction of cts(M T ) [32] ■ However, it is not straightforward to obtain the inclusive 
1=1 component of o (e+e" — > KKn) from the measured value of o (e+e" — > K^K^-k^). In 
fact, using only SU(2)j symmetry this is not possible. There are two further assumptions that 
need to be done in order to obtain the expression usually employed : 



a (e+e" KKn) = 3 a (e+e" -> KsK^) 



(C.5) 



The first one is to assume that the processes e+e" — > KKtx are dominated by i^*-exchange. 
According to recent Dalitz plot analyses, [2H] this is indeed a good approximation. How- 
ever, SU(2)j symmetry and K* dominance do not allow to relate a (e+e" — > KsK~tt + ) and 
a (e+e" — > K + K~tt°) as given by Eq. (10.5[) . Under K* dominance there are two intermediate 
chains that account for each final state : 



A (K + K~tt°) = A (e+e" -> (K*~) K + K'tt K + ) 
+A (e+e" -> (K* + ) K~ ->• K + n° K~) 



B + C . 



A (K-K°n + ) = A (e+e" (K*+) K~ K\ + K~) 

+A (e+e" ->■ (K*°) K° ->■ k + K~ K°^j = V2 (B - C) . (C.6) 
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Accordingly we conclude that the relation a (e + e~ — > K°K~tt + ) = 2 o (e + e~ — > K + K~tt°), 
which is necessary to derive Eq. (IC.5j) . can only be obtained by neglecting the interference 
terms arising from Eqs. (IC.6j) . 

We have checked the accuracy of this assumption using two different parameterizations 
for the involved hadronization amplitudes. First, we have employed a parameterization fol- 
lowing the KS-like model used for tau decays into KKn modes, using the values in Ref. [16] . 
Furthermore, we have used our expressions obtained within R%T in Sect. [3J In both cases, 
we have not set the contributions of other resonances than the K* to zero, although we have 
checked that they are of very little importance. With both kinds of parameterizations either 
at T (r — > KKtcu.,-) or at o (e + e~ — > KKn) the error of assuming that interference effects are 
negligible is at least of order 30 % in e + e~, being even larger in r decays. In the same way, 



Parameterization 


o (e + e~ -)• K°K~n + ) /2 a (e+e~ ->• K+K'ti ) 


m 


0.74 
0.36 



Table C.I: Check of the validity of relation (|C.5j) -and thus of neglecting interference effects- for 
different hadronization parameterizations. The Q 2 -endpoint in the cross-sections has been taken 
at ~ 2 GeV 2 . We estimate the error of the RxT prediction to be around 30%. 



taking for granted the Tightness of the very accurate assumptions of SU(2)i and K* dominance 
in r — > KKitv t , it is still not possible to relate the widths to K + K~n~ and K~K°n° in a 
model independent way. Under the hypothesis of negligible interference one would obtain 

2 T (r K + K-7r'iy T ) = T (r — >■ K~K°tt°i^ t ) . (C.7) 

As can be observed in Table \C.2\ within the above considered models this relation does not 
hold. 



Parameterization 


r (r -> K-K°7T°) /2 T (r K + K--k~) 


m 


0.62 
0.36 



Table C.2: Check of the validity of relation (|C.7j) -and thus of neglecting interference effects- for 
different hadronization parameterizations. We estimate the error of the R^T prediction to be around 
30%. 



D Numerical input 

For the numerics we use, if nothing is specified, the masses given in Ref. [35]. From the 
analyses of Refs. jSTJHT] we use, as input, the following numerical values of the parameters 
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that appear in our study 



F = 0.0924 GeV , F v = 0.180 GeV , F A = 0.149 GeV , 
My = 0.775 GeV , M K * = 0.8953 GeV , M ai = 1.120 GeV. (D.l) 

Then we get A', A" and Ao from the first equalities in Eq. ( 13T?|) . Incidentally we can also 
determine the value of Ma — 0.91 GeV. Notice that this value for Ma is slightly lower than 
the result obtained in Ref. |26j. Our preferred set of values in Eq. ( ID.lj) satisfies reasonably 



well all the short- distance constraints pointed out in Sect. HJ with a deviation from Weinberg 
sum rules of at most 10%, perfectly compatible with deviations due to the single resonance 
approximation. 
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